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Abstract 

We explore the many body physics of a Bose condensed atom gas at finite temperature through 
the Raman transition between two hyperfine levels. Unlike the Bragg scattering where the phonon- 
like nature of the collective excitations has been observed, a different branch of thermal atom 
excitation is found theoretically in the Raman scattering. This excitation is predicted in the 
generalized random phase approximation (GRPA) and has a gapped and parabolic dispersion 
relation. The gap energy results from the exchange interaction and is released during the Raman 
transition. The scattering rate is determined versus the transition frequency uj and the transferred 
momentum q and shows the corresponding resonance around this gap. Nevertheless, the Raman 
scattering process is attenuated by the superfiuid part of the gas. The macroscopic wave function of 
the condensate deforms its shape in order to screen locally the external potential displayed by the 
Raman light beams. This screening is total for a condensed atom transition in order to prevent the 
condensate from incoherent scattering. The experimental observation of this result would explain 
some of the reasons why a superfiuid condensate moves coherently without any friction with its 
surrounding. 

PACS numbers: 03.75.Hh,03.75.Kk,05.30.-d 
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FIG. 1: Feymann diagram illustrating the mediation process: 1) For a plasma two charged exci- 
tations of momentum k and k' mediate their interaction via a plasmon of momentum q; 2) For a 
Bose gas, two excitations with one atom number unit mediate their interaction via a phonon-like 
collective excitation. 

I. INTRODUCTION 

Among the various approximations existing in the literature to describe a diluted Bose 
condensed gas at finite temperature, the generalized random phase approximation (GRPA) 
has been the subject of several studies 

mm 

. This approximation has attracted a 
special attention since it is the only one in the literature with two important properties: 1) 

npriQ 

in agreement with the Hugenholtz-Pines theorem p, 17|, 18|, 19| , it predicts the observed gapless 
and phonon-like excitations; 2) the mass, momentum and energy conservation laws are 
fulfilled in the gas dynamical description. An approximation that satisfies these properties 
is said to be gapless and conserving 

Besides these unique features, the GRPA predicts also other phenomena, namely a sec- 
ond branch of excitations and the dynamical screening of the interaction potential. These 
phenomena appear also in the case of a gas of charged particles or plasma. The possibility of 
a second kind of excitation has been explained quite extensively in . There is a dis- 

tinction between the single particle excitations and the collective excitations. In the case of 
a plasma, the first corresponds to the electrically charged excitations and its dispersion rela- 
tion is obtained from the pole of the one particle Green function. The second corresponds to 
the plasmon which is a chargeless excitation whose the dispersion relation is obtained from 
the pole of the susceptibility function. The plasmon mediates the interaction between two 
charged excitations. More precisely, during the interaction, one charged excitation emits a 
virtual plasmon which is subsequently reabsorbed by another charged excitation (see Fig.l). 
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Remarkably, such a description holds also for a Bose gas with single atom excitations 
carrying one unit of atom number and with gapless collective excitations with no atom 
number. The poles of the Green functions have a similar structure above the critical point. 
But below this critical point, the existence of a macroscopic condensed fraction hybridizes the 
collective and single particle excitations so that the poles of the one particle Green function 
and the susceptibility function mix to form common branches of collective excitations 
Thus, at the difference of a plasma, the presence of a condensed fraction prevents the direct 
observation of the atom-like excitation through the one particle Green function. 

The dynamical screening effect predicted in the GRPA appears much more spectacular 
in a Bose gas. The screening effect of the coulombian interaction is well known to explain 
the dissociation of salt diluted in water into its ions (see Fig. 2a). But it also provides an 
explanation to the superfluidity phenomenon i.e. the possibility of a metastable motion 
without any friction. Most of the literature on superfluidity is usually devoted to the study 
of metastable motion in a toroidal geometry like, for example, an annular region between 
two concentric cylinders possibly in rotation jlol.[lll|. In this simply connected geometry, the 
angular momentum about the axis of the cylinder of the superfiuid is quantized in unit of h. 
The metastability of the motion is explained by the impossibility to go continuously from 
one quantized state to another due to the difficulty to surmount an enormous free-energy 
barrier. This is not the situation we want to address in this paper. We are rather focusing 
on the explanation of the superfiuid ability to fiow without any apparent friction with its 
surrounding. 

The Landau criterion is a necessary but not sufficient condition for superfluidity. It 
tells about the kinematic conditions under which an external object can move relatively 
to a superfluid without damping its relative velocity by emitting a phonon-like collective 
excitation. For a dilute Bose gas at low temperature, it amounts to saying that this relative 
velocity must be lower than the sound velocity [9!]. The external object is assumed to 
be macroscopic and can be an impurity |l2!], an obstacle like a lattice [3] or even the 
normal fluid In particular, this criterion does not taken into account the fact that 
the normal fluid is microscopically composed of thermal excitations. In a Bose condensed 
gas, even though their relative velocity is on average lower than the critical one, many of 
these excitations are very energetic with a relative velocity high enough to allow the phonon 
emission. 
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FIG. 2: Illustration of the screening effect: (a) In water, the interaction force between the ions Na~^ 
and Cl^ of the salt is screened by the presence of water molecules and the coulombian potential 
Vcouii^) is reduced by the relative permittivity factor /C ^ 80. (b) In a Bose condensed gas, a 
similar effect occurs. The condensed and thermal atoms represented in blue and red respectively 
correspond in good approximation to the superfluid and normal fluid. The interaction potential 
V{r) displayed by these thermal atoms on condensed atoms are pictured qualitatively by the green 
line. The macroscopic wave function associated to the condensed atoms deforms its shape in order 
to locally modify the superfluid mean field interaction energy represented by the blue line. The 
net result is a total screening of the interaction potential by this mean field energy, which prevents 
binary collision processes between condensed and thermal atoms. In this way, one can explain 
qualitatively the metastability of a relative motion between the superfluid and the normal fluid. 

In the GRPA where these excitations correspond to the thermal atoms and under the 
condition of the Landau criterion, such a process is forbidden as shown qualitatively from 
Fig. 2b. The effect of an external perturbation of the condensed atoms caused for example 
by the thermal atoms is attenuated by the dynamical screening. This screening is total in 
the sense that no effective mutual binary interaction allows a collision process which would 
be essential for a dissipative relaxation of the superfluid motion. 
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The purpose of this paper is to show that these pecuhar phenomena could in principle 
be observed in a Raman scattering process. This process induces a transition for a given 
frequency uj and a wavevector q determined from the difference of the frequencies and the 
wavevectors of two laser beams jlO]. For each wavevector corresponding to the transferred 
momentum, one can arbitrarily tune the frequency in order to reach the resonance energy 
associated to the excitation. Unlike the Braggscattering which allows the observation of the 



Bogoliubov phonon-like collective excitation 



, |9|, the Raman scattering is more selective. 



Not only the gas is probed with a selected energy transition and transferred momentum, but 
the atoms are scattered into a selected second internal hyperfine level. Through a Zeeman 
splitter, they can be subsequently analyzed separately from unscattered atoms. According 
to the GRPA, the scattered thermal atoms become distinguishable from the unscattered 
ones and thus release the gap energy due to the exchange interaction. In a previous study 
[sl, we showed that this gap appears as a resonance in the frequency spectrum of the atom 
transition rate at q ^ 0. The possibility of momentum transfer allows to analyze the 
influence of the screening of the external perturbation induced by the Raman light beams. 

The paper is divided as follows. In section 2, we review the time-dependant Hartree-Fock 
(TDHF) equations for a spinor condensate and study the linear response function to an 
external potential which gives results equivalent to the GRPA. Sections 3 and 4 are devoted 
to the Bragg and Raman scatterings respectively. Section 5 ends up with the conclusions 
and the perspectives. 



II. TIME-DEPENDANT HARTREE-FOCK APPROXIMATION 

We start from the time- dependant Hartree-Fock equations for describing two component 
spinor Bose gas labeled by a = 1,2. The atoms have a mass m, feel the external 

potential Va6(r,t) and the Hartree and Fock mean field interaction potential characterized 
by the coupling constants gab = 47raafc/m expressed in terms of the scattering lengths ttab 
between components a and b {H = 1). Note that no Fock mean field (or exchange) interaction 
energy appears between condensed atoms. These equations describe the time evolution of a 
set of spinor wave function ?/'a^j(r,t) describing Ni atoms labeled by i and depending on the 



5 



position r and on the time t. For the condensed mode {i = 0), these are: 

+ g - Vu V*, 

For a non condensed mode (z 7^ 0), these are 



-J 




2^ - »'12 
2m 




(2) 



E,(2^?ll|^l,,f + ^?12|^2,,f )iV, ^712 A^.^2,,^l,i 

^?12 A^.A*,,^2,, E, (2^722 |^2,,f + gi2\^l,?)N, 

The non condensed spinors remain orthogonal during their time evolution in the thermo- 
dynamic limit. In general, the spinor associated to the condensed mode does not remain 



orthogonal with the others. But according to 14| , the non orthogonality is not important in 
the thermodynamic limit for smooth external potential. Another way of justifying the non 
orthogonality is to start from an ansatz where the condensed spinor mode is described in 
terms of a coherent state and the non condensed ones in terms of a complete set of orthog- 
onal Fock states i.e. |\t') ~ exp(Ejyo ^i^] ~ '^■'^■) nj^o(4)^' |0) where c| is the atom creation 
operator in the mode i and bj = y/N^Yla I ^^^"^1 j'4'a,o- The theory remains conserving 
because the conservation laws are preserved on average but becomes non number conserving 
since the quantum state is not an eigenstate of the total particle number operator. This 
procedure is justified in the thermodynamic limit since the total particle number fluctua- 
tions are relatively small during the time evolution. In contrast, instead of using spinor 
wavefunctions, the alternative method based on the use of excitation operators is number 
conserving 3|, |j] . 

The atom number Ni for each mode is supposed time-independent in the TDHF. Strictly 
speaking, a collision term must be added in order to allow population transfers between the 
various modes. These equations are valid in the collisionless regime i.e. on a time scale 
shorter than the average time between two collisions r ~ 1 / {aahnvx) where vt = f3m is 
the average velocity and aab = Svra^j, is the scattering cross section. In these conditions, the 
resulting frequency spectrum has a resolution limited by Aa; ~ 1/r. The magnitude order 



of resolution of interest is given by the gab^^s so we require Aoo /gatn ~ ^/O'lb'n/ Pgatn <^ 1 
which is generally the case when a^^n <^ 1. These conditions are fulfilled for the parameter 
values considered in this work. 

In the following, we will restrict our analysis to a bulk gas embedded in a volume V. At 
t < 0, we assume all atoms in thermodynamic equilibrium in the level 1 and that Vat = 
except for V22 = ^0 which is constant and fixes the energy shift between the two sub-levels. 
In that case, the solutions of the TDHF are orthogonal plane waves with i corresponding to 
the momentum k: 

' -AS _ "''Pl'tk-r - <ifi)l ( 1 \ 

i-Sl)^ ^ [0) 

where we define the Hartree-Fock energy for atoms with momentum k: 

efk = Ck + S-iilSn - no(5k,o) (4) 

where ek = k^/2m and where the condensed and total particle densities are no = Nq/V and 
n = ^k/V^- Eq-dlD corresponds to the dispersion relation of the single particle excitation. 
At equilibrium, 

K = A^k(l - 5k,o) = l/(exp[/3(efk^ - /.)] - 1) (5) 

is the Bose-Einstein distribution. Below the condensation point, the chemical potential 
becomes /i = cq = gii{2n — uq) and the macroscopic occupation Nq is fixed to satisfy the 
total number conservation. 

For t > 0, we apply an external potential. For the Bragg and Raman scatterings, these 
are respectively: 

Vii = Vb cos(q.r — ujt) (6) 
^12 = V^ijexp[i(q.r-cjt)] (7) 

We solve the system through a perturbative expansion: 
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The equations of motion for the first order corrections are for the case of Bragg and Raman 
scatterings respectively: 



2m 



^2 



Vu +$^^7ii(2 - 5k',o5k,o)(V^S'V^;i + c.c.)iVk, 



■41 



(9) 



(10) 



These two set of integral equations can be solved exactly using the methods developed 
in Defining the Fourier transforms: 



/ Jo 



one obtains in the level 1 for the condensed mode: 



for the non condensed modes (k 7^ 0): 



(11) 



(12) 



■f 2m /C(q', uj'){uj' + iO - et+q' + et) 



(13) 



and in the level 2 for all modes: 

duj' 



/c 
■( 



27ri 



■X 



(14) 



/Ci2(q', u;')(^' + «0 - c^o - Ck+q' + Ck + (25'ii - 5'i2)n + (5k,o(fi'i2 - fi'ii)^o) 
These formulae resemble the one obtained from the non interacting Bose gas excepted for 
the mean field term in f[T^ and the extra factors representing the screening effect. For the 



Bragg scattering, these factors can be written as 



3|: 



A(q,cu) 



tu + iOf 



/C(q,^) 



A(q,cu) 



{uj + iOy - Cq + 25'ii"'oeq 



where 



^(q,'^) = (1 - 25'iiXo(q,tu))[('x^ + ^0) -eq ] - 85(iiXo(q,tu)fl'ii?^oeq 
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(15) 
(16) 

(17) 




FIG. 3: Diagrammatic representation of the scattering of an atom by an external potential 

is the propagator for the collective excitations, = i^c^q^ + e'j^ is the Bogoliubov excitation 
energy, c = ^/gunojm is the sound velocity and 

1/ ^ u; + 20 + ek - ek+q 
is the susceptibility function describing the normal atoms. For the Raman scattering, it is 

^i2(q, uj) = l- 5'i2Xo,i2(q, t^) (19) 

where 

Xo,i2(q, ^) = -r; — —7\ . ^— ^ — -I ^ — (20) 

V/ ^ + «0 - cuo + Ck - Ck+q + (25(11 - gi2)n + dk,o(5'i2 - 911)^0 

Knowing the Fourier transform of the potential Vii,q/,a;' = X]± ^^B^q',±q/2(co'' + uj) and 
^i2,q',i^' = '^VR5q,q//(u;' + iO — uj) , Eqs. fll2lll3|[T^ are calculated using the contour integration 
method over uj' by analytic continuation in the lower half plane. As a consequence, the poles 
of the integrand tell about the excitation frequencies induced by the external perturbation. 
The pole of the propagator containing k corresponds to atom excitation involving one mode 
only while the poles coming from the screening factors correspond to the excitations involving 
all modes k collectively. Thus, the TDHF approach predicts both single atom and collective 
excitations. Note that the single mode excitation is not possible for the condensed atoms 
since the corresponding pole is compensated by a zero coming from the screening factor. 
The expressions ( |T2]|13|I14I) have an interpretation shown in Fig. 3. An atom of momentum k 
is scattered into a state of momentum k + q' by means of an external interaction mediated 
by a virtual collective excitation of momentum q'. 
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III. BRAGG SCATTERING 



Let us first review the Bragg scattering process . Up to the second order in the Bragg 
potential, the atoms number for any mode k can be decomposed into an unscattered part: 



j^unscat ^ 



1+ / dh{i;f^i;fl + c.c.] 
Jv 



(21) 



and a scattered part: 

AT--^* = I d'v\ip\'^l\' (22) 



'V 

Instead of evaluating the second order term, A^^"«c<^* is determined through the conservation 
relation = ]Sf^^^'^°-^ -)- A^^™*. Generally speaking within the sublevel 1, the scattered 
atoms cannot be distinguished from the unscattered ones. But in order to understand 
the underlying physics, we assume that distinction is possible. Within the second order 
perturbation theory, the quantity of interest is the scattered atom rate per unit of time and 
is expected to reach a stationary value after a certain transition time. In the following, we 
shall analyze these transition rates for time long enough that transient effects disappear. In 
these conditions, a perturbative approach is still valid for very large time provided that the 
scattered atom number remains low compared to unscattered ones. This last requirement 
is always satisfied with a sufficiently weak external perturbation. 

At zero temperature, only the condensed wave function is modified and Eq. f|T2l) becomes 
after contour integration over uj': 



V^S!i(r,t) = f^S(r,t)$:e±- 



(23) 



The response function is only resonant at the Bogoliubov energy ie^. Also no transient 
response appears at zero temperature. Using ( |T2l) and ( |T3i) . the total number of scattered 
atom can be obtained by determining the total momentum: 

P = E /" dh k— ^i,k = E / dh l^fipq (24) 
^ Jv I k Jv 

In the large time limit, the total momentum rate is related to the imaginary part of the 
susceptibility response function x = x' ~ "^x" through 

dP t^oo _ / Vb ^ 



W-Ht^ /(q,^) (25) 
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Using Eq. (l23|) . we recover that: 

X"(q, oo) = 7rS^No{6{LO - ej) - 6{iO + ej)) (26) 

where = Cq/eq is the static structure factor. The delta function comes from the relation 
6{x) = lim^^oo sin(xt)/(7rx). The result (!26l) obtained in the GRPA is identical to the one 
obtained from the Bogoliubov approach where can be calculated equivalently from the 
four points correlation function p, [l^ . But in any case the generated phonon like excitation 
is still a part of the macroscopic wave function tpi^o{r,t). 

At temperatures different from zero, the poles become imaginary which means that any 
Bogoliubov excitation is absorbed by a thermal atom excitation 1^. This phenomenon 
is known as the Landau damping. So for long time, only the residues of f|T2l) with poles 
touching the real axis contribute whereas the others give rise to transient terms negligible 
for long time. Thus the perturbative part becomes: 

*S;i(r. t) |^*S(r. t) J: \ (27) 

(1) t oo Vb (0) ( e^*'^* g-i(ek±q-ek)t \ giiq.r 

= ^4>,^, t) X: (^(^^ - ^(^q^ ^^^^ _ j _ ^^^^ ^ ^^, (28) 

Using the property A(q, u) = A*(— q, —ui), the total number in the condensed mode reaches 
a constant value 

A,--J»fM^M±^ (29) 



V 2 y |A(q,cj^l2 
and the scattered thermal atom rate is given by: 



rft V 2 y ^ |/C(q,cj 

From (l25!l . we deduce for the imaginary susceptibility: 

X"(q,^) = -— Imf-p-^) (31) 

The basic interpretation of these formulae is the following. At finite temperature, the col- 
lective excitation modes created by the external perturbation are damped over a time given 
by the inverse of the Landau damping. So the number of collectively excited condensed 
atom reaches the constant value fl29p when the produced collective excitations rate compen- 
sates their absorption rate by thermal atoms. This constant value is higher for a transition 
frequency and a transferred momentum close to the resonance a; = ec ~ ie'^ 



q ■ 
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The formula (l30!l is a generalization of the Fermi- Golden rule when the screening effect 
is taken into account. The external potential perturbs the thermal atoms of momentum k 
in two channels by transferring a momentum ±q and a transition energy ±uj such that the 
resulting single atom excitation has a momentum k ± q and a kinetic energy e^-i-q = ± cu. 
The presence of the screening factor amplifies or reduces the scattering rate. Amplification 
(or ant i- screening) occurs for a frequency close to the resonance energy ec of the collective 
excitations. On the contrary, dynamical screening occurs for a frequency close to the pole of 
the screening factor and is total for transition involving condensed atom aX u = Cq. Thus, 
in GRPA, attempt to generate incoherence through single condensed atom scattering is 
forbidden at finite temperature. Only collective excitations affect the condensed mode but 
they are damped and therefore cannot contribute to effectively transfer condensed atoms to a 
different mode [sl]. It is taught in standard textbooks that, in the impulse approximation 
used for large q, the response of the system is sensitive to the momentum distribution of the 
gas, since the atoms behave like independent particle. In particular, a delta peak is expected 
to account for the presence of a condensate fraction. The difficulty of the observation of 
this peak could be explained by this impossibility of a single condensed atom excitation at 
finite temperature. For completeness, let us mention that interaction with thermal atoms 
can be also totally screened and inspection of the formulae (fT6|) shows that this happens for 
eg = i^Cq — c^q^ j^. Fig. 4 shows these features in the frequency spectrum for the total 
momentum rate (|3Tll at fixed q. We choose the typical density observed experimentally for 



87 



Rb at the trap center 



These results can be put in direct relation with the analysis of impurity scattering 15|. 
Indeed, the dynamic response function is related to the dynamic structure factor through 
the fiuctuation-dissipation theorem: ^(q, u;) = x"(q, Ci;)/7r(l — exp(— /^cj)). The dynamic 
structure factor is directly connected to the transition probability rate P(q, cu) that an 
external particle or impurity changes its initial momentum p and energy Ep into p + q and 
E'p+q = Ep + uj respectively: 

V{q,uj)=2iT\V^\^S{q,uj) (32) 

where Vq is the Fourier transform of the interaction potential between the impurity and the 
atom gas. The total rate of scattering Fp results from a virtual process involving emission 
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FIG. 4: Imaginary susceptibility x" of a bulk Bose condensed gas for eq = 27r x 30kIIz versus the 
detuning frequency 5uj. The superfluid fraction is 94%, gun = 2it x 4.3kHz, kBT/gun = 2.11 
and afiU = 5.6 10~^. The black dashed/solid curve is the rate calculated in absence/presence of 
the screening factor. Both regimes of screening and anti-screening are displayed close to the zero 
eg and to the resonance Cc respectively. In particular, the screening prevents the observation of a 
huge delta peak associated to the condensed mode. 

and absorption of the collective excitations: 

Tp = J2^7r\V^\'Siq,E^+^-E^) (33) 
q 

= E^^l rr ^''^(^i^ + - ^^^^ - ^k-q)K(l + K-J (34) 

q,k P+'l ~ 

As a consequence, the impurity scattering is possible provided that the energy and momen- 
tum are conserved in a effective collision with a thermal atom of momentum k mediated by 
a virtual collective excitation. Note that total screening prevents impurity scattering involv- 
ing ongoing and outgoing condensed atoms. In contrast, for temperature close to zero, the 
Landau damping approaches zero since xoi^.,^) — ^ so that the application of Eq. (l2^ to 
(!33|) leads to an on-energy shell process of absorption and emission of a collective excitation. 
We obtain: 

Tp = J] 27r|VqP^q« + 5±,+ )5i±4 + i?p+q " £^p) (35) 

±,q 



where = l/(exp(/3eq) — 1). This Umit case leads to the apparent interpretation of an 
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impurity interactin g w ith a thermal bath of phonon-hke quasi-particle. This situation has 
been considered in [ig] in the study of the impurity dynamics. Instead, Eq. flMl) provides a 
generahzation for higher temperature emphasizing that any external particle can excite a 
single thermal atom alone but not a condensed one. 

IV. RAMAN SCATTERING 

The conclusions so far obtained in the Bragg process can be extended straightforwardly 
to the case of Raman scattering with the difference that only one channel of scattering is 
possible. For the purpose of simplicity, we choose the case g = gab- Also this channel is 
easier to access experimentally. Defining the detuning 6u = u — ujq, exphcit calculations of 
the spinor component (fMl) in the second sublevel give: 

'^'^ ' V^i2(q,^^) /Ci2(q,cJo + ek+q-ek-5(n)y i(5c^ + ek + 5(n-ek+q) 

So we obtain for the atom number in the mode k: 

c?A^2,k t^oo _ , .2^('^^ - ek+q + Ck + gn)N]^ 

—— = 2'kVr — — — (37) 

dt |/Ci2(q,cj)|^ 

By summing over all the modes, we obtain the density rate transferred in level 2: 

^ = 1(E ^-^^1^) ^^nX'uici, uj) (38) 

k 

where we define the imaginary part Xi2 = Xu " ^Xi2 of the intercomponent susceptibility 
function: 

Xi2(q, ^) = Xo,i2(q, ^)/ (1 - ^Xo,i2(q, ^)) (39) 

n 

This last formulae is also the one obtained in the GRPA [4]. Again we find a similar structure 
as the intracomponent case. In this process, thermal atoms with an initial momentum k 
and energy e^^ = 2gn + ek are transferred into a second level with momentum k + q and 
energy ^2,k+^ = 9^ + ^k+q provided 6u! = ^2,k+(i ~ ■ absence of screening, a resonance 
appears at the detuning = eq — gn. The first term corresponds to the usual recoil energy 
while the second is the gap energy gn that results from the exchange interaction. During 
the Raman transition, the transferred atoms become distinguishable from the others and 
release this gap energy. 
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FIG. 5: Imaginary susceptibility Xi2 of a bulk Bose condensed gas eq = 27r x 30Hz versus the 
detuning frequency buj. Parameter values are the ones of Fig. 4. Left and right graphs represent 
the same curves but the right graph is in logarithm scale. The black dashed/solid curve is calcu- 
lated in absence/presence of the screening factor while the dotted curve represents the Bogoliubov 
approximation. See the grey curve for a magnification of the black solid curve (x25) 

The scattering rate is determined through the imaginary part of the susceptibility Eq. (15^ 
versus the transition frequency uj and at fixed q. Figs. [5] and [6] show the corresponding 
resonance around this gap in absence of screening. The screening effect strongly reduces 
the Raman scattering and, in particular, forbids it for atoms with momentum k such that 
k.q = 0. This case corresponds to buj = eg and includes also the condensed atoms (k = 0). 
The graphs illustrate well the effect of the macroscopic wave function that deforms its shape 
in order to attenuate locally the external potential displayed by the Raman light beams 
and to prevent incoherent scattering of the condensed atom. The experimental observation 
of this result would explain some of the reasons for which a superfluid condensate moves 
coherently without any friction with its surrounding. Anti-screening occurs in the region 
close to the resonance frequency ec of the collective mode. At zero temperature, we recover 



= eq [l7l] while for non zero temperature the collective modes become damped for q 7^ 
These results can be compared to the one obtained from the Bogoliubov non conserving 

n 

approximation developed in p and valid only for a weakly depleted condensate. This 
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FIG. 6: Idem as Figl5]but for eq = 2tt x 300Hz. Here, the broadening of the curves is much more 
important. 

approach implicitly assumes that the only elementary excitations are the collective ones 
and form a basis of quantum orthogonal states that describe the thermal part of the gas. 
Consequently, this formalism predicts no gap and no screening. Instead, the intercomponent 
susceptibility describes transitions involving the two collective excitation modes of phonon 
and of rotation in spinor space ek: 

no , 1 V- ulJn^ + 5±,_) 



45: 



— Cq + iO V ^ c<j + iO ± e 



(40) 



±,k ■ ^ ^k±q 

where w-t^k = iK^k + fi'^o)/2e^ ± 1/2]^/^. This function does not preserve the f-sum rule 
associated to the SU{2) symmetry. In contrast to the GRPA, a delta peak describes a 
spinor rotation transition of the condensed fraction, and two other transitions involve the 
excitation transfer from a phonon mode into a rotation mode and the excitation creation 
in the two modes simultaneously. For small q, these processes remain dispersive since the 
frequency transition depends on the momentum k. As a consequence, the resulting spectrum 
shown in Figs. 5 and 6 is broader. In particular, the process of creation in the two modes 
favors transition with positive frequency. Note also the maximum of the curve separating 
the region involving a transition atom-atom like (high k) and the one involving a transition 
phonon-atom like (low k). 

All these features established so far for the bulk case allow a clear comparison between the 
GRPA and the Bogoliubov approaches. In the real case of a parabolic trap, the inhomogene- 
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ity induces a supplementary broadening of the spectrum that prevents the direct observation 
of the screening. This effect as well as the finite time resolution and the difference between 
the scattering lengths will be discussed in a subsequent work. 

V. CONCLUSIONS AND PERSPECTIVES 

We have analyzed the many body properties that can be extracted from the Raman 
scattering in the framework the GRPA. The calculated spectrum allows to show the existence 
of a second branch of excitation but also the screening effect which prevents the excitation 
of the condensed mode alone. 

The observation of phenomena like the gap and the dynamical screening could have 
significant repercussions on our microscopic understanding of a finite temperature Bose 
condensed gas and its superfluidity mechanism. On the contrary, the non-observation of 
these phenomena would imply that the gapless and conserving GRPA is not valid. In 
that different approximation has to be developed in order to explain what will be 

observed. As an alternative, the idea to use the Bogoliubov approach has been also discussed. 
But unfortunately, the violation of the f-sum rule is a serious concern regarding this non 
conserving approach [5]. All these aspects emphasize the importance of the experimental 
study of the Raman scattering at finite temperature. 
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and from the German AvH foundation. 
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